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^ ' We discuss price variations distributions in foreign exchange markets, characterizing 

them both in calendar and business time frameworks. The price dynamics is found to 

be the result of two distinct processes, a multi- variance diffusion and an error process. 

The presence of the latter, which dominates at short time scales, leads to indeterminacy 

C^ , principle in finance. Furthermore, dynamics does not allow for a scheme based on 

■^T ■ independent probability distributions, since volatility exhibits a strong correlation even 

• ' at the shortest time scales. 
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' O ' 1. Introduction and definition of time 

Q ' The characterization of price dynamics in financial markets is an old but still 

puzzling problem. At the beginning of the century, Bachelicr W proposed to con- 
sider price variations as independent realizations of identically Gaussian distributed 
variables, while in the '60 Mandelbrot introduced symmetric Levy stable distri- 
butions. In the last years, Mantegna and Stanley (0], and see also Rpl) provided 
t;:j- I evidence that a Levy stable process well reproduces the central part of high fre- 

CO . quency price variation distribution, while the tails are approximately exponential. 

The Bachelier's Gaussian shape is recovered only on longer time scales which are, 
O^ ' typically, of order one month. The common point of ||l|-[£| is that price dynamics 

0^ ■ is considered as the result of independent random variables. This kind of approach 

seems inadequate, since there is evidence of volatility correlations on long run [|6|- 

@- 

Indeed, price distribution strongly depends on how one measures time flow. The 
'"O ' choice of time index is twofold, calendar time and business time. Business time is 

^ ■ the sequence of integers n = 1, 2, 3, which indexes successive established quotes. 

fj , These integers correspond respectively to the calendar times ^1,^2, Therefore, 

calendar time (which is monotonically increasing) is a stochastic process of business 
time. The relation can be inverted by considering n as a function of calendar time, 
rN ' i.e. n — nt, but in this case the function is defined only on the sequence of quoting 

j^ ' calendar times. The price dynamics, therefore, can be described with respect to 

business time by means of the series S'„ or with respect to the calendar time by 
means of St = Sn^ ■ In the first cases quotes S are defined on all integers and lags 
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are all equal, in the second, they are defined only on quoting times and lags are 
unequal. 

The main result of this paper is that a price quote is due to two distinct inde- 
pendent stochastic processes: an error process superposed to an underlying price 
process. The latter evolves following calendar time, while the former is due to an 
erroneous evaluation of a market operator, and its natural frequency is marked 
by business time. The resulting price variation distribution is, therefore, the con- 
volution of two distributions associated to these two distinct processes: the error 
process distribution, which does not change with time, and the distribution of the 
underlying process, which scales, on the contrary, with calendar time. 

The error process produces always a gap between two consecutive price estima- 
tions even if they are almost contemporary. This phenomenology strongly reminds 
quantum mechanics, where a measurement result always has a minimum uncer- 
tainty as stated by the Heisenberg Principle. Following this comparison, one can 
state an indeterminacy principle in finance: a price is never given with a precision 
less than a natural constant for that market. 

Another important fact which is due to this phenomenology is that two consec- 
utive price variations cannot be considered fully independent random variables, but 
they exhibit a very peculiar anticorrelation as we will see. Moreover, we provide 
evidence that volatility is so strongly correlated that remains substantially constant 
inside the largest lag we consider, and therefore the usual Levy stable scheme seems 
to be not appropriate. 

In this paper we examine the high frequency price variation distribution of three 
foreign exchange markets, the Deutsch Mark/US Dollar (DEM-USD) exchange in 
1993 (1,472,140 quotes) and in 1998 (1,620,843 quotes), the Japanese Yen/US Dollar 
(JPY-USD) in 1993 (570,713 quotes), and the Japanese Yen/Deutsch Mark (JPY- 
DEM) in 1993 (158,878 quotes). The quotes represent the value of one US Dollar in 
Deutsch Marks and Japanese Yens in, respectively, the DEM-USD and JPY-USD 
cases, while they represent the value of one Deutsch Mark in Japanese Yens in the 
DEM-JPY case. The price changes are given in pips, which indicate a DEM/10, 000 
in the DEM-USD case, and a Yen/ 100 in the JPY-USD and JPY-DEM cases. All 
the data sets analyzed in this work have been provided by Olsen & Associates. 

We deal only with bid price, since we have found more spurious data in the ask 
price distributions due to wrong transcriptions. This is probably related to the fact 
that in recording process quotes are specified by bid price and the last two digits 
of the bid/ask spread. Nevertheless, all the results can be fully reproduced if one 
deals with ask price or average price. 

The paper is organized as follows. In Section 2 it is shown that a price variation 
is given by the composition of two distinct stochastic processes and that the hy- 
pothesis of independent consecutive price variations is incorrect, due to the presence 
of both price anticorrelation and volatility correlation. In Section 3 the price vari- 
ation distribution for a single business lag and for the minimum calendar lag (two 
seconds) are computed. The latter turns out to be, basically, the error distribution. 
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In Section 4 we provide evidence that the underlying price variation distribution at 
a generic calendar time lag is given by a symmetric Gaussian process whose stan- 
dard deviation is itself a Log-normal process. In Section 5 sonic final remarks are 
reported. 

2. The indeterminacy principle in finance 

The aim of this Section is to show that the set of reported bid quotes is a real- 
ization of stochastic process which is the composition of two processes whose origin 
and meaning is very different. The first is the ordinary multiplicative stochastic 
process which determines the evolution of underlying price, while the second is a 
superposed stochastic noise which somehow accounts for the erroneous evaluation 
of the underlying price by the market operators. 

In order to demonstrate our point it is better to consider the price evolution in 
business time. 

Let us assume that it exists an unobserved underlying price S'„ which accounts 
for the real relative value of two currencies and which follows the ordinary evolution 
rule 

Sn+l — Sn + Rn (2-1) 

where, as usual, the i?„ are, at the lowest approximation, identically distributed 
variables which: a) have vanishing average, b) only depend on the stochastic calen- 
dar lag Ai„ = tn+i — tn, c) are uncorrelated, i.e. the average of the product RnRm 
vanishes for n y^ m, d) their typical size (standard deviation) is proportional to Sn- 
The last property it is necessary to ensure that the process is multiplicative on the 
large time scale. Notice that we do not assume that the i?„ are independent, and 
the reason for that will be clear at the end of this Section. 

Let us also assume that the observed price (which is the recorded quote) slightly 
differs from the underlying price because of an erroneous evaluation of the operator. 
The relation between the underlying and the observed price will be 

Sn = Sn + En (2.2) 

where the En are zero mean, uncorrelated identically distributed variables. Inde- 
pendence, again, is not assumed. At variance with the i?„ they are independent 
from calendar lags Ai„: in fact they are a consequence of price evaluation process, 
and therefore they arc present at each business time. 

As a result of the joint action of this two processes one has that after m business 
lags the observed price evolves according to 

m—l 

Sn+m = Sn + En+ni ^ En + / ^ Rn+i (2.o) 

i=0 

Now one can appreciate that price variation differently depends on two con- 
tributions, the first, which represents the evolution of the underlying price, is the 
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Fig. 1. DEM-USD 1998 exchange rate: variance (2.4) (circles) compared wi th a hnear fit 2A+Bm 
with A = 5.9 ± 0.2 and B = 1.85 ± 0.01, neighbouring autocorrelation (2.5) (slan ting crosses) 
compared with —A from the previous linear fit, non-overlapping au toco rrelation (2.6) (rhombuses) 
compared with zero, and square root of quadratic autocorrelation ( p.7| ) (crosses) compared with a 
scaling fit Cm in the range 20 < m < 50 with C = 2.60 ± 0.03. 



sum of the m variables, the second, which represents the uncertainty inherent to 
the quoting operation, is always the difference of only two uncorrelated identically 
distributed variables. 

If this theoretical framework is valid there should be various consequences. First, 
according to Eqn. (p^) the variance of a price change after m business time should 
be 



<(^, 



n'\-m 



Snf >= 2A + Bm 



(2.4) 



where < • > indicates the average over the probability distribution, and A —< E^^ > 
and B =< R^ >. 

This last average also runs on all possible calendar lags between two successive 
quotations. The occurrence of this property can be appreciate in Fig. 1 which refers 
to DEM-USD 1998 exchange quotes, and price variations are given in DEM/10, 000. 
Results are highly consistent with the linear behaviour with A — 5.9 ± 0.2 and 
B = 1.85 ±0.01. 

The non-vanishing value of A has an important meaning, since it suggests that 
market price models based on continuous time approach [|l4| are not appropriate: in 
fact, this would at least require in this case a vanishing price change in the limit of 
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vanishing lag. In this case, price change variance should be zero in the limit n ^ 
(i.e. A = 0), while in Fig. 1 it can be clearly appreciate the non occurrence of this 
fact. 



The second, more peculiar, consequence of Eqn. (2.3) is that the neighbouring 



autocorrelation of two consecutive price variations after m business time is 

<^ (*5r!+m ~ '-'n)('-'n ~ '-'Ti-m) >— ^A (^-5) 



where A is the same parameter of (2^). The similar autocorrelation for non- 



overlapping price variations {I > 0), according to Eqn. (2.3), are 



< {Sn+l+m — Sn+l)iSn — Sn-m) >— (2.6) 

Indeed, these last two equations represent a strong test in order to verify the good- 
ness of the proposed price modelization, since phenomenology is very peculiar and 
crudely differs, for example, from Markov anticorrelations. 

Both these relations are well satisfied as it can be appreciated in Fig. 1. The 
same results can be found by using other data sets: A = 7.7±0.2 and B — 2.59±0.01 
for DEM-USD 1993 exchange rate, A = 8.4 ±0.4 and B = 2.96 ±0.03 for JPY-USD 
1993, and A = 1.2 ± 0.2 and B = 4.88 ± 0.01 for JPY-DEM 1993. In the last two 
cases price variations are given in Yen/100 (pips). 

It should be noticed that in all cases 2 A is about ten times larger then B. This 
means that the observed price changes are largely dominated by the error effect on a 
short time scale. Only after about ten business lags, the underlying price variation 
is comparable with error and it can be partially appreciated. Ten business lags 
correspond, approximatively, to two minutes in calendar time for the DEM-USD 
cases. This means that price change distributions on time scales of few minutes are 
deeply affected by error, and do not reflect the real price variation. 

At this point we are sufficiently convinced that a price is a fuzzy variable, not 
only because of the obvious bid/ask spread, but especially because the bid price 
itself cannot be given at each time with an absolute precision less than ^/A. We 
call this simple fact the indeterminacy principle of markets. 

In conclusion of this Section, let us consider that Eqn. (^.6|), which states that 
the autocorrelation of two non-overlapping price variations is zero, may lead to the 
idea that these variables are independent. The inconsistence of this hypothesis can 
be provided by computing the following non-overlapping quadratic autocorrelation: 

< {Sn+m+l — Sn+l) {Sn ~ Sn-mj > — < (Sn — Sn-m) > (2-7) 

If two non-overlapping absolute price variations were independent, one simply would 
have that this quantity vanishes. 



On the contrary, the quadratic autocorrelation (2.7) is sensibly different with 
respect to zero, as it can be appreciated in Fig. 1, where its square root is plotted. 
In fact, at large business times m, it shows a linear behaviour Cm. The numerical 
estimations give C = 2.60 ±0.03 for DEM-USD 1998 exchange rate, C = 3.07 ±0.08 
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for DEM-USD 1993, C = 2.7 ± 0.1 for JPY-USD 1993 and C = 4.2 ± 0.3 for JPY- 
DEM 1993. 

The last result gives an important information concerning the quadratic auto- 
correlation of underlying price process: 

< R„R„^^ > - < R„ > ~C (2.8) 

independently on m at least for business time lag m in the range 20 < to < 50 (about 
ten minutes for DEM-USD 1998 exchange rate). In other words, Eqn. (2.8) asserts 
that absolute price have an autocorrelation which is substantially independent on 
the time separation to.. Indeed, this is a very surprising result, and it confirms the 
presence of a strong correlation on absolute price changes not only on time scales 
from minutes [Esl to months |E2|,E3|, but also on shortest time scales of seconds. 



Indeed, absolute price exhibit a very reach phenomenology also for what concerns 



long term behaviour. In particular, recent results on long memory |12|l^ show 
power-law correlations up to time scale of one year for |A5|^, where jAS*! is price 
variation and a; is a real number, whose exponent depends on x (multiscaling). 

3. Single lag distribution and minimum calendar lag distribution 



According to (2.3) the process changes are determined by the joint effect of 



underlying evolution and added noise. The equation (2.3) can be easily rewritten 
in calendar time as 

St+T = St + En{t+T) ~ En{t) + Rt,T (3.9) 

where t and t + t are two calendar times in which quotes are established, and Rt.r 
represents the variation of underlying price from time t to time t-\-T corresponding 



to the sum in (2.3) 



In consequence of (3.9) and taking into account the uncorrelation of all the 



variables, the probability distribution of a price change R is given by the convolution 

P^{R) = Q{AE) (g) Pr{R) (3.10) 

where Q{AE) is the error distribution and Pt{R) is the probability distribution of 
underlying price R. It should be noticed that Q{AE) is directly the probability of 
the differences of the two variables En{t+T) ~ £'«(*) and not of a single one. The 
probability Q{AE) does not depend on the calendar time lag r, at variance with 

PriR). 

The error distribution Q{AE) and the underlying distribution Pt{R) cannot be 
directly observed, since the observable probability distribution is only the Pr{R). 
Nevertheless, the analysis of the price change variance ( p.4D suggests that at a single 
business time the distribution is largely ruled by the error distribution. Further- 
more, since two seconds is the minimal calendar lag, (i„+i — 1„ is always a multiple of 
two seconds) the corresponding price variations are a subset of those corresponding 
to the minimal business lag n ~ \. At the light of these informations, we expected 
that the minimum calendar time lag distribution Pt is practically equal to Q. 
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Fig. 2. DEM-USD 1998 exchange rate: price change distribution at minimum calendar lag of 2 
seconds (circles) and at 1 business lag (crosses) in log-linear scales as functions of price variation 
measured in pips (DEM/10, 000). The minimum calendar lag distribution basically coincides with 
the error distribution tanh(7/2) exp(— 7I A_E|), and this gives the parametric fit 7 = 0.40 ± 0.01 in 
the business time range —23 < m < +23. 



In Fig. 2 we plot in log-linear scale the two seconds calendar lag distribution. 
The evidence of Fig. 2 is that the probability distribution Q{AE) is exponential of 
the form 

(3.11) 



Q{AE)^ tanh- exp(-7|A£;| 



for six orders of magnitude, and the numerical fit of the parameter gives 7 = 
0.40 ± 0.01 in the business time range —23 < m < -f 23. 

A simple check of self-consistence can be made by computing the variance of 
the error distribution Q{AE) and comparing the result with 2A of Eqn. (2.4). In 
fact, Q{AE) is the probability distribution of the difference of two independent 



variables E, 



n(t+T) 



-E„(f), each one having variance equal to A. After some algebra, 



the variance of distribution ( 3.11 ) reads 



< (AE)^ > = 



1 



cosh 7 — 1 



(3.12) 



By inserting into this formula the estimated 7 = 0.40 ± 0.02, one has < (AE)'^ >- 
12 ± 1, in agreement with the value we have found in the previous Section 2 A - 
11.8 ±0.4. 
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In fig. 2 the minimum business lag distribution is also plotted, for comparison 
with the two seconds distribution. The central parts of the two distributions are 
substantially the same for three orders of magnitude, in this region the error process 
rules. Outside the region the single business lag distribution exhibits more persis- 
tent tails. These large events, corresponding to calendar time lags larger than two 
seconds, are due to underlying price variations instead of erroneous price changes. 

4. Distributions at different calendar time lags 

In this Section we discuss the probability distributions of price variations for 
calendar lags of length r from its minimum of two seconds up to time scales of 
order of few minutes. In particular, we try to retrieve informations about the 



underlying price distribution PriR) in Eqn. (3.1C). 

A possible strategy is to guess a functional parametric expression for Pr{R), and 
then to compare its (5(Ai5)-convolution (3.10|) with our experimental probability 



distributions. In order to make a reasonable hypothesis on Pr{R), let us come back 
to business time approach and write down the underlying price change i?„ as the 
product of a volatility an and a Normal Gaussian variable with zero mean and 
unitary variance W„: 

Rn = (Jn Wn (4.13) 

Notice that this decomposition always allows for independent Normal variables Wn ■ 
From Eqn. (p^) and independence of Wn , one can derive a similar relation for the 
volatility: 

< ^,V'+™ > - < ^l >'^ C' (4.14) 

for m up to 50. 

The simplest interpretation of the last result is that volatility remains sub- 
stantially constant in business time lags of length m, but it assumes different val- 
ues for temporally far lags (otherwise the last expression should vanish). Than 
C^ ~< (jfj^ > — < cr^ >^. In other words, in a business lag of m steps the volatil- 
ity process is substantially frozen, while one has m independent realizations of the 
Gaussian process. This result identically holds for calendar time, where one says 
that in a time lag of length t up to, about, ten minutes for DEM-USD 1998 ex- 
change rate, one has a constant volatility. This volatility can obviously change with 
respect to that of another time lag of length r temporally far. 

In other terms, volatility itself is a stochastic process with a characteristic time 
much larger than about ten minutes. The underlying price variation Rt^r on a 
calendar lag of length r is then the sum of several independent Gaussian variables 
with same variance. Therefore, it can be written as: 

Rlr = (JtV^Wt (4.15) 

where at is constant in the calendar lag (t, t + t) and Wt is a Normal Gaussian 
random variable. 
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The price process framework is completed when an exphcit expression for the 
volatiUty probabiUty distribution is given. With this aim, let us recall recent results 
PpJ^] , where it was found that returns probability distributions of a stock market 
daily index and of a foreign daily exchange are symmetric Gaussian distributions 
whose standard deviations are themselves Log-normal processes. Borrowing this 
result and taking into account the previous considerations, the following expression 
for Pt{R) comes out: 

PriR) = / da L{a) / dr G„^[r) (4.16) 

where L(cr) is a Log- normal distributions of parameters /i and lo 

-{^\\a~ [if 



L{a) = y= exp 



2uj^ 



(4.17) 



and G^r^{r) is a Gaussian distribution with zero mean and standard deviation 
a/tit, integrated on the unitary interval centered on the integer R. 



We have performed a numerical computation of Eqn. (4.16) via a Monte- 
Carlo approach, and then we have convoluted the result with the error distribution 
Q{/S.E). For the DEM-USD 1998 exchange rate we have found a parameters choice 
(/i = —1.5 and to — 0.8) which gives a good agreement between the guessed and 
the experimental price change distributions, as shown in Fig. 3, for three values of 
time lag r: 4, 40 ± 2 and 400 ± 20 seconds. In order to have a reasonable number of 
occurrences, the experimental distributions are computed with a 5% tolerance on 
r. 

The consistency of the hypothesis of a Log-normal volatility which is constant 
during a whole time lag can be directly checked computing ^< af^ > — < a"^ >'^ 
for such a distribution and comparing the result with the experimental value of C . 
In the DEM-USD 1998 case one has that the first computation gives ~ 2.77 against 
the value C = 2.60 ± 0.03. Therefore the constancy of volatility seems reasonably 
confirmed. 

The conclusion is that the underlying volatility itself evolves following a stochas- 
tic process, with characteristic times larger than about ten minutes, while on shorter 
time scales the snapshot variance is mainly due to the error process. The hypoth- 
esis of a Log-normal underlying volatility which seems once more confirmed by 
experimental data, could be the result of a multiplicative process. 

5. Conclusions 

We have found evidence that two distinct stochastic processes, an underlying 
process and an error process, are present in price change dynamics in foreign ex- 
change markets. The latter process, which can be put in connection with the price 
estimations of market operators, suggests a quantum like nature for price changes 
variables, in the sense that they have an irreducible intrinsic indeterminacy. For 
this reason market models based on continuous time limit seem to be not adequate. 
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Fig. 3. DEM-USD 1998 exchange rate: price change distributions for calendar time lag r equal to 
4 seconds (circles), 40 ± 2 seconds (crosses) and 400 it 20 seconds (rhombuses) in log-linear scales 
as functions of price variation measured in pips (DEM /IO. 000). The lines repr esent the numerical 
computation of the guessed distribution Pt{R) (1.16) convoluted with Q{AE) (3.11) for fi = —1.5 
and LO = 0.8. In the inset the small price changes range is magnified for the cases 4 seconds and 
400 seconds. 



Another main consequence of the indeterminacy principle is that price changes 
cannot be considered the result of independent stochastic processes, because of 
the presence of correlations even at the shortest time scales. In fact, the spurious 
anticorrelations in the observed price changes |lq ] can be easily explained in the 
light of the error process. 

At this point it seem quite clear that a moving chart over few observed prices 
is able to give a more accurate estimations of the underlying price at a given time. 
Nevertheless, the number of observed prices in the moving chart cannot be too large 
(i.e. calendar time lag cannot be too long) otherwise too large price variations due 
to the underlying process would decrease accuracy. 

We expect to find effects of the indeterminacy principle also in high frequency 
market quotes of single stocks. On the contrary, the error component in change 
distributions of composite market indexes are probably not so important because 
of the averaging of errors due to the large number of different stocks involved. 

Finally, the volatility of the underlying process exhibits a strong autocorrelation, 
since it turns out to be substantially constant up to time scales of, at least, ten 
minutes. This is a clear evidence that the underlying volatility is a stochastic 
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process with a larger characteristic time, while variance at short time scales is 
basically due to the error process. We have proposed that the underlying volatility 
follows of a Log-normal process and that the underlying price change distribution 
is a symmetric Gaussian whose standard deviation is the Log-normal stochastic 
volatility. This hypothesis is in a very good agreement with experimental data. 
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